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ABSTRACT 


This thesis will examine the equations of motion for a 
spinning ballistic projectile. The goal of such an 
examination is to determine the possible mechanisms by which 
a directed energy weapon may induce sufficient instability as 
to significantly alter the projectile's flight path. A 
ballistic projectile is generally launched with a "fire and 
forget" philosophy. The desired impact point is determined 
before firing. It may be possible to alter the projectile in 
such a way that it fails to follow the desired trajectory 
thereby missing the intended target. Several variables appear 
to be worthy of investigation to assess their contribution to 
a required instability or range reduction. Skin friction drag 
may be increased from surface roughness generated by a pulsed 
energy source. The results that this thesis will examine 
include: impulse generated by the laser interaction, 
additional Magnus effects and aerodynamic drag. Moment 
induced instability may also result from these in the form of 
a Magnus moment or drag torque. Increasing the drag force 
appears to be the most promising theoretical solution to 
defeating an incoming spinning ballistic projectile. 
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TABLE OF VALUES 


c s 



specific heat of steel 

0.15 cal/gr-C 

d 



155mm projectile diameter 

.155 m 

g 



acceleration due to gravity 

9.8 m/sec 2 




spin decelerating coefficient 

.08 

k d 



drag force coefficient 

.004 

k h 



damping moment coefficient 

1.2 

k l 



lift force coefficient 

1.0 

% 



over turning moment coeff. 

2.5 

k n 



normal force coefficient 

1.3 

k t 



Magnus moment coefficient 

0.3 

L 



155mm projectile length 

70 cm 

l m 



latent heat of melting 

65 cal/gr 

m 



155mm projectile mass 

45 kg 

t m 



melting temperature of steel 

1450 C 

u 



wind speed over the projectile 

600 m/s 

u s 



speed of sound in air 

350 m/s 

V 



speed of the projectile 

600 m/s 

CD 



155mm projectile spin frequency 

220 hz 

Pair 



standard density of air 

3 

1.2 gr /cm 

P steel 



density of steel 

3 

8 gr /cm 

* Many of these terms are approximations based on diagrams and 
conversations with aerodynamics researchers at the Army Research 
Laboratory in Aberdeen, Maryland. The exactness of these values is 
not significant to the magnitude estimations in this thesis. 
Additional special case values will be introduced in the text. 
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I. 


INTRODUCTION 


This thesis is directed toward solving a stated goal of the 
Defense Planning Guide. Our current armed forces are vulnerable to 
gunfire attack due to an antiquated counterbattery system of 
equipment and technique. Specifically, in today's battlefield, 
surgical precision is required and may involve protecting civilian 
lives in a confined suburban environment. There is a need to 
develop a new system using state of the art technology to destroy 
incoming projectiles minimizing loss of life. Laser technology has 
given us the ability to track, classify and engage incoming 
ballistic rounds to meet this goal. This thesis will explore the 
use of laser energy to accomplish the defeat of the projectile. 

This mission responds to the Department of Defense "Defense 
Planning Guidance (DPG) for FY 1996-2001" dated 23 May 1994. This 
involves countering armor, mortar and mobile artillery threats to 
U.S. forces conducting an amphibious invasion, U.S. forces 
defending a coastal region from invasion and to coastal cities and 
safe zones under U.S. protection during peace keeping operations. 
There are several specific guidelines in the DPG: 

- "U.S. and coalition forces would have several key objectives in a peace 
enforcement or smaller - scale operation, each of which would require certain 
capabilities:... Establishing and defending zones in which civilians are 
protected from external attacks." 

- "The first priority in defending against a large-scale attack will most 
often be to neutralize the enemy's offensive capabilities." 

"U.S. military capabilities will be configured to achieve required 
objectives while minimizing American casualties." 
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A. DEFEATING THE LAUNCHER BEFORE PROJECTILE IS A THREAT 


There exists a need for a Counter-Battery Weapons System 
(CBWS) to neutralize the enemy's ability to launch large caliber 
ballistic projectiles at U.S. forces and civilians under U.S. 
protection. The CBWS needs to be able to determine the launchers 
position from one in-coming round and surgically destroy that round 
in as rapid a manner as possible. A secondary purpose that this 
thesis will not address is eliminating the launcher based on cuing 
information from the primary system. 

The Combat Systems Science and Technology Curriculum at the 
Naval Postgraduate School has a course sequence designed to solve 
these kinds of problems and propose system solutions. The class of 
June 1995, examined the problem solution of defeating projectile by 
preventing the launch. The class considered this problem to be 
"shooting the archer" as opposed to "shooting the arrow." A 
complicated and expensive proposal was developed using acoustic and 
electromagnetic sensors to detect an initial firing and then a 
guided projectile to destroy the launcher and nearby personnel 

B. DEFEATING THE INCOMING PROJECTILE 

This weapon system should have an integrated, multi-unit and 
autonomous capability to designate the target and automatically 
engage the threat using a laser tracking beam and a higher power 
laser destruction beam. Countering the threat from armor, mortars 
and mobile artillery, the system must be able to respond against 
projectiles ingressing at up to three km/s from a range of up to 
120 km. The response time for incoming projectile detection 
through battle force threat warning and counterbattery action must 
therefore be less than 45 seconds. The system should be capable of 
distinguishing friendly fire from enemy fire through command 
functions and will also be able to project the most probable impact 
point of the in-coming round to determine if impact could cause 
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major damage to friendly forces and to provide warning of imminent 
danger to ground forces. Projectiles that pose no threat will not 
be engaged. The initial goal of the CBWS will be to achieve a hard 
kill of the projectile. 

Modern weapon R&D is producing long range over-the-horizon 
(OTH), high caliber ballistic weapons systems which are small 
enough to be highly mobile. These weapons systems are being 
deployed in areas which are becoming an important part of the U.S. 
military strategy in which Naval forces are being used to support 
ground forces. The shift from emphasis on open ocean warfighting 
to the littoral environment as outlined in "Forward...From the Sea" 
will increasingly place forces at sea within range of these weapons 
systems. Working with the Marine Corps and their Operational 
Maneuver From the Sea (OMFTS) warfighting doctrine will require 
support in various roles. These include: suppression of enemy 
gunfire, neutralization/denial of enemy operating areas, Naval 
Gunfire Support (NGFS) , interdiction/neutralization of reinforcing 
elements, and evacuation protection. These capabilities will 
require execution in day or night, all weather, and OTH scenarios 
to counter the threat. Specifically, the threat must be engagable 
in an area lOOnm wide and 7 5nm deep (from an offshore range of 
25nm) . This must be accomplished with an accuracy of several 
meters to minimize collateral damage and fratricide. The purpose 
of the laser weapon system is to provide the capability to achieve 
these objectives. 

Within these requirements, our current inventory of naval 
gunfire support (NGFS) weapons is inadequate to counter this 
threat. The technological gains made in other Navy systems has 
occurred while naval guns have remained virtually unchanged. The 
maximum range of our most capable system, the 5" 54-caliber M 45, 
falls far short of the desired ranges and accuracies. 
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C. ALTERNATIVE SYSTEM SOLUTIONS TO DEFEAT THE LAUNCHER 


There are few alternatives to developing this system. One 
non-material alternative may be to adapt existing counterbattery 
systems (Firefinder radar) for use by naval forces. Because the 
precision required to surgically neutralize enemy weapon systems 
does not presently exist in the Navy's arsenal, improved tactics 
cannot lead to a solution of the stated problem. 

Currently, both the Army and the Marine Corps have ballistic 
tracking systems in service to locate hostile guns. These again, 
fall significantly short of the required accuracy and do not 
address the problem of locating and destroying the incoming 
projectile. Nevertheless, these systems providing a search 
capability used in conjunction with a more accurate fire 
control/targeting system may provide good initial or secondary 
queuing prior to the terminal phase of the desired weapons system. 

Another alternative involves improvements to the Mk 45 
projectiles. Several options are being explored including 
extending the ranges of the projectiles with advanced conventional 
solid propellants such as liquid propellant or electrothermal 
chemical propellant. Further improvements to the projectile could 
provide a "smart" seeker which would aid in the terminal portion of 
the ballistic trajectory. Use of the existing TOMAHAWK cruise 
missile as a weapons delivery system would require modifying the 
missile to be able to target a point defined by latitude and 
longitude vice using predetermined scene maps and an image matching 
method. But, again, these systems would do nothing to protect the 
forces and civilians from the projectile itself. 

D. THE SOLUTION APPROACH OF THIS THESIS 

The problem can then be summarized as "how it is possible to 
defeat an incoming ballistic round?" The flowchart below describes 
the major solution areas and subareas. This thesis is directed 
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toward solving the problem of deflection/destruction in flight. 
Specifically, the question "what can be done to change the 
ballistic projectile's dynamic trajectory?" This is accomplished 
by examining the forces that control the dynamics and what forces 
may be applied to change the dynamics (Figure 1). 



Single Pulse Multiple Pulses Drag yaw Torque Magnus Effects 


Figure 1 Defining the problem is the first step in determining 
possible solutions to defeating an incoming ballistic projectile. 

The solution of providing armor to defeat the round is 
considered elementary and could be considered as limiting the scope 
of the other solutions to larger caliber rounds such as the 
155mm. The ability to defeat an incoming round in flight, and 
ultimately the success of the proposed system, is explored through 
several theoretical calculations involving laser interaction with 
the surface of a projectile. Two main areas of research appeared 
promising in the problem solution: ablative impulse and induced 
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surface roughness. The impulse problem was handily dismissed due 
to the large energy requirements. However, the induced surface 
roughness may have a significant effect on a ballistic projectile's 
stability and range. This instability desired to deflect the round 
is most easily obtained by altering the Magnus moment of the 
projectile. This instability, however, proved to be of 
insufficient magnitude to defend the projectile's target. The 
additional drag created by the increased surface roughness appears 
to be a promising solution. 

This thesis estimates the magnitude of effects that could de¬ 
stabilize or reduce the range of a projectile in flight. Most 
previous work on exterior ballistics has focused on the means for 
achieving stable flight and maximum aerodynamic efficiency. The 
magnitude of the effects we seek may demand great pointing accuracy 
and other laser CBWS problems may arise. These are not addressed 
here. 
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The 155mm projectile is used as a model to conduct these 
calculations because it meets the requirements of a medium to long 
range artillery threat and is a common weapon in many international 
forces. The 155mm projectile has many variations of casing design 
and explosive warheads that are used. The most common projectile 
is the 155mm HE M107 which is used as a general purpose round by 
the United States and most other countries that fire the 155mm 
ballistic projectile. The HE M107 round weighs 95 pounds and is 
27.5 inches long. The casing is composed of a medium carbon 4340 
steel (a carbon, nickel, monel alloy) that varies in thickness from 
.25 to .75 inches at the nose. 



Figure 3 An exaggerated description of the increasing angle of 
the projectile with respect to the trajectory along the 
trajectory. 
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As the projectile travels through the air, it spins at a 
frequency of 220-250 hz. The spin is imparted by rifling within 
the bore of the launcher. A variation in the spin rate is due to 
the differing charges that may be used to fire the projectile. A 
typical exit velocity of 600 meters/second will be used throughout 
the calculations. In addition to spin, the ballistic projectile 
also undergoes a coning motion about the axis of flight with an 
origin near the center of mass(18.1 inches from the nose tip). The 
angle of this coning motion is typically 3-5 degrees as the 
projectile is reaching its apex and then increasing somewhat as the 
arc of trajectory curves downward toward the target. The coning 
rate is between 5 and 20 hz. 

Ultimately, the system can only be successful if the 
theoretical results can be further refined and tested 
experimentally to verify expected system performance. The purpose 
of this writing is to find possible theoretical means of using 
laser energy in a CBWS to achieve the national goals stated in the 
DPG. Specific details of system design will not be addressed. 

The thesis will present arguments for five possible mechanisms 
of achieving a "kill" of the projectile. All of these "kills" will 
be achieved by deflecting vice destroying the projectile. The 
first theoretical method discussed to defeat the projectile will be 
momentum transfer by directed energy impulse. A more complex 
argument discusses the stability equation for a spinning projectile 
and addresses the most likely component, the Magnus effect, in 
relation to achieving instability. Finally, increased aerodynamic 
drag is examined as the most likely means to achieve a successful 
engagement against an incoming spinning ballistic projectile. 
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II. IMPULSE DELIVERED BY LASER INTERACTION 


A. LASER EQUATIONS AND IMPULSE THEORY 


Before examining any calculations of ablative impulse, it is 
necessary to introduce the equations which will be used. The 
equations that will be used were developed in class notes used at 
the Naval Postgraduate School for various courses examining 
directed energy capabilities (Schriempf, 1974) . The terms that 
will be used are defined below: 


F Q = (l - SI) I The laser power density absorbed where 31 is the 
Reflectance, I is the laser power density arriving 
at the surface and F Q is the power transmitted 
through the surface per unit area. 


t 


P 


The pulse length (duration). 


D The thermal diffusion length is the distance 

hr 

required for the temperature to drop l/e of the 
surface value. 


P 


T 


0 


m 


c s , Cl 
L m' L v 


The material density. 

The initial temperature. 

The melting temperature. 

The vaporization temperature. 

The specific heat of the solid, liquid. 
The latent heat of melting, vaporization. 
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The boiling temperature. 


K „ 

K = — C 

P 


The thermal diffusivity is equal to the product of 
the thermal conductivity and the specific heat 
divided by the density. 


The criterion for vaporization of a diffusion thickness is 
given by: 


F o ~n ~ P ^s( T m~ T o> +L m +c l( T b~ T m> + V ■ 


( 2 . 1 ) 


When numerical values are inserted, the L y term dominates the 
equation simplifying to the form: 

T, 




( 2 . 2 ) 


Since, D d « 2 /TkO , this then becomes: 


F O — 2 -*!= 1 V P 


(2.3) 


Solving this equation for the pulse length gives 

2 

■^v 2 

T p >4K — P . 


(2.4) 


Now to solve the specific impulse on a target will require 
several assumptions before this can be carried any further. 

1) The beam radius, R, is much larger than the diffusion 
length. 

2) The target is very thick compared to the diffusion length. 

3) The melting will be ignored. 
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The time for the surface to reach vaporization temperature is: 


x b = 




(2.5) 


The surface erosion speed is taken from equation 2.1 
substituting A z for the diffusion length. This represents the 
thickness of the material evaporated during a time At. The 
erosion speed is then expressed as: 

'-y 1 /[C s( T m -T 0 )*C p (T v -T m )*I. v ) . <2 ' 6) 


Since L m <L y ; Cj«C s =C, then U — 1/ [ L y +CT V ] . 


(2.7) 


The vapor blow-off speed is derived from the momentum 
conservation, p v U v = p s U s , resulting in: 


p v U v = F Q / (L V +CT V ) . 


(2.8) 


The associated pressure with the blow-off is then: 

2 R /ft 

P = P V U V = p v — (gas constant/molecular weight) T y so, U v = 

and combining this with equation 2.8 results in: 

P = F 0^]f^ / {L V +CT V ) * F 0 JCjT^/[j3(L v+ CT v )] (2.9) 

since C«3— for most metals. 

M 

Now, the specific impulse delivered during an energy pulse is: 
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J in p ^ x p x b^ • 


(2.10) 


Substituting equations 2.5 and 2.9 into 2.10 and using typical 
values results in a simplified form: 

X V , 2, (2.11) 

I m = 8£ 0 (l-6-|), [(taps) dyn-s/cm ] 


where E Q is the pulse energy per unit area (fluence) 
= F 0 t p , T p units are |is. 

This equation will generally result in agreement with 
experimental data within a factor of two. This is sufficient for 
the estimates of this thesis. If the results of further 
calculations indicate a possible effect within a factor of two from 
current technology, then better estimates would be in order. 

One additional aspect of laser theory should be addressed 
here. If the E Q is very large, the delivered impulse will fall off 
with increasing values of E Q . This is due to the absorption of 
laser energy within the vapor near the surface of the target. At 
some value of power density, the plume of vaporized material blow- 
off from the target will become hot enough that the surrounding air 
or vapor are sufficiently ionized to absorb large amounts of the 
laser radiation. This phenomenon is called laser absorption. 

Another problem is called beam decoupling. When the laser 
beam electron density increases to a value where the plasma 
frequency is larger than the laser frequency, then the light 
interaction with the free electrons dominates. The laser is cut¬ 
off from the surface. 

With these in mind, it is important to find the critical 
energy density at which these processes occur. The upper limit 
density is about 10 23 electrons per cubic cm electron density where 
the material reaches a temperature of full ionization. But, this 
is not accurate enough since the ionization degree is very 
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sensitive to the temperature. Figure 4 shows the predictions for 
the impulse as a function of the pulse energy (Nielsen, 1972) . 
These values will be used later in the calculations for ablative 
impulse. 



Figure 4 Specific impulse delivered to solid targets by 1.06 
micron laser radiation (Nielsen, 1972). 
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B. ABLATIVE IMPULSE CALCULATIONS 

l. Single Pulse Impulse 

The impulse calculations are examined first since this would 
be the easiest and most direct method of deflecting an incoming 
round. A Laser induced ablation from a single pulse would be the 
preferred method of destruction if it were physically possible. 
This would require a less sophisticated system that could "fire and 
forget." The problem is complicated when a system is required to 
accurately track and reengage a target multiple times. This 
requirement also detracts from the ability to rapidly engage the 
next subsequent target in a dense battlefield. With this in mind, 
and using the equations for critical energy of a laser system, it 
is possible to determine the feasibility of the single pulse 

system. 



14 




To further simplify the initial calculation, a perpendicular 
impulse will be assumed. If the required energy is at physically 
attainable levels, then further calculations for reduced angle 
impulse will be in order. 

The impulse transferred by ablation is found by equation 2.12 
to be proportional to E Q for Tp^lO|is and E Q >20J/cm when the term 
in the bracket becomes small compared to one: 

I m = 8 E o< 1 -( 6 ^f >) • 

E 0 ( 2 . 12 ) 


The impulse is calculated from the change in momentum required 
where the momentum of the moving projectile is, 


p=mv=1.8x10* gem/s. 


(2.13) 


For a single pulse system, a shipboard close-in defense would 
require a deflection on the order of 50 meters for every 1000 
meters traveled by the projectile. The result would be a .05 rad 
change in direction or a change in momentum of 


(Ap) re gd = * 05 ^ ■ 8xl ° 9 ) =l0 8 gcin/s. 


(2.14) 


Then, the impulse transferred in a 10 square cm area to the 
projectile is estimated by: 


dm* regd : 


(Ap ) regd i 


=10 dyns/cm‘ 


(2.15) 


Finally, referring to equation 2.12, the energy fluence required 
will be 10 °J/cm with a pulse duration of 10 6 seconds Presently, 
this is not within the realm of possibilities for current laser 
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systems. Therefore, further calculations of reduced angle impulse 
are not required. However, there is another way to examine the 
problem. It may not be a close-in defense system that is engaging 
the target projectile. Without placing any limitations on the 
amount of deflection required for a given length of travel, a 
calculation may be conducted for a change of yaw angle. A 
significant change of yaw angle of about 10-20 degrees would change 
the drag significantly and alter the range of the projectile. 
Using a spinning top approach and a small angle change, a first 
approximation again rules out the possibility of sufficient energy 
for a single pulse. 

2 

From Figure 4, values of Energy 30 J/cm and I/E of 5 dyn s/J 

2 

show that an impulse of 150 dyn s/cm could be attained. 

To determine the required impulse, one finds the required 
torque from the required change in angular momentum in the 
following steps: 


the torque, T=Pa, 


(2.16) 


1 2 

the moment of inertia, I = —MR , (2.17) 

and the angular momentum, L = Jco = 10>d? (2.18) 



Figure 6a The torque applied. Figure 6b A change in angle. 





Therefore, the torque can be expressed as: 


dL 

dt 


= 10) 


dn 

dt 


= 1 0) 


da 

dt 


F-a . 


( 2 . 20 ) 


The impulse per unit area required is: 


^Preg 


F-At 

A 


JO) 

IFa' 



JO) 

aA 


•Aa. 


( 2 . 21 ) 


Using the values from Figure 4, E czit = —with the result of 30 
MW. This is an energy fluence far beyond the laser absorption 
threshold and is an implausible solution with current technology. 

2. Multiple Pulse Impulse 

This leaves the less desirable option of impulse delivered by 
multiple pulses to be examined. Using the results of equation 
2.21, the total energy required must be divided into optimal energy 
pulses. If E opt = 30—^, then 10 5 pulses are required. Using 10 ps 
pulses at optimal energy, it is possible to determine how much 
engagement time is possible as the projectile cones. The 
engagement time is estimated by the following using a coning rate 
Of 5 Hz: 


U 

G) c = 2 —Hz also v=R c G) c = atana G) c . 


( 2 . 22 ) 
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With the established coning rate and angular velocity, an 
engagement can only be achieved within a period of .06 seconds. 
Applying these pulses over the engagement time for a 1 ns pulse 
results in a power requirement of 6xl0 10 Watts which is still well 
outside of an expected capability. 

The impulse calculations demonstrate that it is currently not 
possible to achieve a "kill" on a ballistic projectile with the 
ablation energy transfer alone. The next step is to determine how 
the stability may be affected by changing the aerodynamic forces 
acting on the projectile. 
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m. STABILITY 


There are several topics that should be explored before 
beginning the calculations to prove or disprove that a viable CBWS 
can cause instability. Most importantly, a sound understanding of 
spinning projectile stability theory is established as a framework. 
Then, using the stability as a guideline, one term (the Magnus 
moment) in the final stability equation appears to be most 
promising for further analysis in causing instability. The Magnus 
force and ultimately Magnus moment must be defined and discussed to 
determine whether these components are worthy of further 
calculations to support an instability argument. 

The conditions for stability of a spinning ballistic 
projectile are given by McShane et al. The following discussion is 
a summary of the results of this analysis with the details given in 
Appendix A. 

A. THE EQUATIONS OF MOTION 

The focus of this analysis is the derivation of the equation 
of motion of the yawing motion due to the various aerodynamic 
forces acting on the projectile. The velocity vector, 3, of the 
projectile with respect to the ground can be decomposed into a 
component along the projectile axis of symmetry, u-j_, and a complex 
component, i; = u 2 + iu 3 , where u 2 and u 3 are the velocity components 
perpendicular to the projectile axis in a coordinate system fixed 

to the projectile. The tangential velocity of the center of mass 

F 

is U. The ratio, is the "vector of yaw" which has the 

magnitude of sin5 . A new variable is defined as: 
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p 

i fvdp 

X.Ae ° , '3-D 



where P is the arclength along the trajectory in calibers and v is 
the spin rate per caliber of travel. The yawing equation of motion 
derived in Appendix A is a second order inhomogeneous differential 
equation for the quantity, X. 

The solutions to the homogeneous equation are given as 
exponential functions in P as: 


X = 


M 


i/i ) ± <w*' aj K- > ' 0+ivA(1±0) /BJ dp 


(3.2) 


where the integrand is a function of the relevant aerodynamic force 
coefficients. The yawing motion is then given by the superposition 
of the two solutions of the homogeneous equation plus a particular 
solution of the inhomogeneous equation: 

X = + L 2^2 + *p’ (3-3) 

The coefficients L ± , L 2 would be determined by the initial yaw and 
rate of yaw on the trajectory. 

The projectile is stable along the trajectory if the solutions 
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k 1 and X 2 decrease in magnitude so that in the limit of long 
trajectory, only the particular solution remains. This means that 
A is independent of the initial conditions. The particular 
solution need not be small however. For large spin, the rate of 
precession and the direction of its axis changes only slowly. As 
the trajectory curves, the yaw for a stable projectile can still be 
significant. In fact, the projectile becomes more stable along the 
trajectory as the velocity of the projectile slows at a greater 
rate than the spin. 

In reality, the yaw will be small if the projectile is stable 
and the solution of the particular integral is small; this is 
almost always the case unless high-angle fire is used, then the 
particular solution may not be small. 

The analysis of stability of the yaw motion concentrates then 
on the two solutions A 1 and A 2 which can be written in the form: 

r i 

f [ — (a 1 +a 4 ——) ±r] dp 


where a 1 , a 4 , r are functions of the aerodynamic coefficients 
(defined in Appendix A) and the spin rate which is not time 
independent along the trajectory. 

One can see that the yaw will decrease if and only if the real 
part of the first bracket under the integral is negative and in 
absolute value greater than the real part of the second. 

In the following, we will summarize this analysis with the 
thought in mind to determine which of the aerodynamic coefficients 
must be changed by laser action, for instance, and by what 
magnitude in order to destabilize the projectile. 


21 



B. THE STABILITY EQUATIONS 


The purpose of this section can be achieved now by defining 
the terms of stability which are derived more thoroughly in the 
appendix. The stability factor S and the stability terms are 
usually expressed as: 

1 k 


Hid ^ 

S 1 = J N +k -2 J H~ J D~ J A~Y~ ' 


ind ^ 

S ^2 = 2Jfj~2J]j-2Jj ,———, and 


(3.5) 


S 


3 


2k 2 J h +(2J t - 2J a ) 


md 2 

A 


The conditions for 
inequalities: 


stability are expressed as the 


1 

S 


< 


S 2 S 3 


and > 0 . 


(3.6) 


It is worth noting that the spin term, v, appears only on one side 
of the inequality. This is useful for determining the effect of 
varying spin rate on a projectile's stability. 

These inequalities yield the desired framework for comments on 
projectile stability. Namely, there are conditions which do not 
satisfy the inequalities that are cause for instability. When 
either of the terms of equation 3.6 become negative, the conditions 
for stability are violated. These are the conditions for 
exploitation in the calculations. It becomes clear in analysis 
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(see the discussion at the end of Appendix A) that the only term 
which can be reasonable exploited to cause instability without 
catastrophic damage to the entire projectile is the Magnus moment 
term, J T . This is true because this term has a polarizing effect 
on the first condition of 3.6. Since this term contains opposite 
signs in each term of the numerator, then it will have an effect in 
opposite directions of the over-all stability. Further, because 
the Magnus coefficient can either be positive or negative depending 
on its relation to the center of mass location, it can change signs 
within each term (but, change them both in opposite directions 
simultaneously). 

The drag term also appears capable of causing a change in the 
stability term, S 2 , since it is of the same order of magnitude as 
the Magnus term. However, the J N term contains the drag 

coefficient as well. So a change in the drag will result in a 
corresponding change in J N' nullifying any effect it would have on 
the stability. In examining equations 3.5, it can be seen that 
increasing the drag will cancel out a similar increase in 
Applying any change in the drag term to either the denominator or 
the numerator of equation 3.6 will result in similar cancelling of 
the effect. This leaves the Magnus term as the only term that can 
be altered in such a way to violate the stability criteria 
expressed in equation 3.6. 

C. THE MAGNUS FORCE AND MAGNUS MOMENT 

The Magnus moment is the most likely component of stability to 
be exploited. By increasing the Magnus effect sufficiently, one 
can cause instability. Before making the argument for instability, 
it is appropriate at this point to discuss the Magnus force and 
Magnus moment. 
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First, the coordinate system of the projectile should be 
established. The vector components of this system will be used to 
define the Magnus terms. A ground-based, right-handed coordinate 
system will be used with a center of mass of the rigid projectile 
serving as the origin. 



The projectile is considered a solid body of revolution. 
Along the axis of rotational symmetry, a unit vector X is chosen 
positive in the direction toward the nose from the center of mass. 
Therefore the total angular momentum of the body can be expressed 
as the sum of: 

1) the angular momentum about X 
and the 

2) total angular momentum about an axis perpendicular to X. 
Both vectors are in the ground-fixed coordinate system. 

The sum of all the forces and moments will determine the 
position of the projectile. For this discussion, however, only the 
notation that is necessary to understand the Magnus contributions 
will be introduced including several additional variables that 
provide an understanding of this effect. 

For the angular momentum about X, the magnitude of the product 
AN, where A is the moment of inertia of the body about X, and N is 
the axial spin (angular velocity with a positive direction of a 
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right-handed screw about x expressed in radians per second) result 
in the total angular momentum represented by the vector AN5c. 

An aerodynamic force or moment (torque) is produced by the 
interaction of the rigid body and the fluid of the atmosphere 
around it. There are several in addition to the Magnus components 
that are considered necessary to form a complete model of 
symmetrically spinning projectile behavior. These are outlined in 
Appendix A. The drag force will be introduced in more detail in a 
later chapter as a model for examining the effect of a roughened 
projectile skin on an aerodynamic coefficient. 

The Magnus effects on spinning projectiles were first 
documented by G. Magnus in the 1800's. He observed that a spinning 
musket shot experienced drift in the direction of spin as if there 
were an additional sideforce acting on the projectile. It has been 
consistently recognized that this effect is an important 
contribution in determining stability and trajectory of spinning 
projectiles, missiles and bombs. Although there has been a great 
deal of interest in these effects, there are few studies that 
examine these effects in great detail. Because of the complexity 
of this effect we will at this stage use the basic principle of the 
Magnus force and Magnus moment as outlined below make a simple 
estimate of the effect of varying surface roughness. To date, 
there have been no models developed to describe the complexities of 
an asymmetric flow caused by an asymmetric body, although one paper 
addresses the non-linear effects of the Magnus force for large yaw 
(Platou, 1956) . 

It is widely accepted that the Magnus effect is produced 
entirely by viscous effects. An asymmetric boundary layer (around 
a symmetric rigid body) is produced by both the spin of the 
projectile and its angle of attack. Most of the research to date 
has examined bodies with low angles of attack and very smooth 
surfaces. A simple illustration shows the asymmetric boundary 
layer and the resulting forces for the low angle of attack. 


25 




In Figure 10, we consider a cross section through the 
projectile with its plane spanned by the trajectory and the 
direction given by £ * X, where € is the vector along the 
trajectory. The bottom view represents the projectile cross 
section in this plane with a boundary layer and the cross velocity 
component of the wind as it adds or subtracts to the boundary layer 
velocity due to the projectile spin in this plane. 

On the upper portion of the cross section, the fluid flow over 
the projectile is parallel to the rotational motion, and on the 
lower side they are in opposite directions. On the upper side, the 
mean velocity in the boundary layer is larger than on the lower 
side. Therefore, the pressure is larger on the lower side and the 
effect is a net pressure force acting upward in the Figure. The 
projectile is acted upon by a net force to the left when viewed 
along the trajectory. 

For a large angle of attack (large yaw in excess of 20 
degrees), the boundary layer begins to separate and vortices flow 
into the outer inviscid flow. The result is an equal and opposite 
circulation causing an additional sideforce. Currently three 
dimensional models of boundary layer behavior are too difficult for 
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Figure 11 Vortices shed due to the Magnus effect and may lead 
to additional force component in the Magnus terms (Power, 1974). 

theoretical predictions beyond the crudest approximations. 

Since the effects on the 155mm projectile are calculable 
within a low angle of attack, the further analysis of the three 
dimensional boundary layer are not conducted in this work. Such 
circulation and vortex shedding may contribute to rapid growth of 
instability as precession occurs and should be further investigated 
in subsequent work. 

For the symmetric body before a laser has induced surface 
roughness, the Magnus force can be fairly accurately described. 
The Magnus force is the result of the fluid interaction of the air 
stream flowing over the boundary layer of a yawed spinning body. 
The magnitude of the force is pd 3 JCpNvsin6, where K p is the 
dimensionless coefficient of the Magnus force. The force is 
perpendicular to the yaw plane. The Magnus force is represented by 
the vector for positive values of N. The result is illustrated 
in Figure 12. 




Some research has been conducted to determine the value of the 
coefficient and relates it to the length to diameter (L/D) ratio of 
the body for low angle of attack and very smooth surfaces in 
supersonic flight (Power, 1974). Since we are trying to determine 
the effect of varying the roughness, this solution is not helpful 
in the end, but may be used to determine a ballpark value of the 
unperturbed coefficient. The results of this research are 
reproduced in the figure on the following page. 
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Figure 13 Angle of attack vs Magnus Force coefficient 
(Power, 1974). 

Although it is the Magnus force that results from the 
roughness that we hope will ultimately cause instability, it is the 
moment arm placement of this force on the body that causes the 
torque about the center of mass. This is the contribution of the 
Magnus moment. If the Magnus force does not pass directly through 
the center of mass of the projectile, then a resulting Magnus 
moment is produced. The magnitude of the moment is given as is 
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the dimensionless coefficient of the moment. This coefficient is 
approximated in McShane, et al. with a dependence upon the center 
of mass location in the projectile (the location of the CM is on 
the horizontal axis of Figure 14) (McShane, 1953). 



(McShane, 1954). 


The Figure depicts firing data obtained for a 4.5 inch rocket 
and although not precise for our purposes, it provides a general 
idea of how the magnitude of the coefficients varies. The value of 
the coefficient is read from the graph at the point where the 
moments are applied in relation to the center of mass position). 
For the 155mm projectile, the center of mass would be near the 
number two on the axis of Figure 14. This gives a general idea of 







representative values for the Magnus moment if a Magnus force were 
applied to the forward third of the projectile (e.g. a Magnus force 
applied at position 4.5 would yield a Magnus moment coefficient 
near zero). 

The value of the moment coefficient may be positive or 
negative as the mean position of the force moves away from or 
toward the center of mass. The hope is to achieve a roughness that 
will produce a force with a different length moment arm from the 
center of mass and thus produce a change in the Magnus moment 
coefficient. The actual center of mass will not have changed, but 
the position of the Magnus force relative to the center of mass 
will have the effect of changing the center of mass for predicting 
the Magnus moment coefficient with the above Figure. The Magnus 
moment lies in the yaw plane and is perpendicular to x. The vector 
[x* (x*v) ] has magnitude vsin5 . The moment is depicted and 
expressed in Figure 15. 

It is hoped that inducing a Magnus moment by increasing 
surface roughness would become a technologically feasible 
possibility. If it becomes possible to use surgical precision (to 
within centimeters) to produce a roughness at a desired location 
from the center of mass, then the resulting torque would be changed 
by the amount r X F, where r is the distance from the center of 
mass and F is the magnitude of the Magnus force. Ultimately, then 
the problem of defeating an incoming projectile could be achieved 
with a calculated amount of energy directed toward a specific point 
on the projectile with predictable results. This would be the most 
efficient use of the proposed system. 
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MAGNUS MOMENT 



1 

MAGNUS MOMENT = M T = pd K t N[xX(xXv)] 


Figure 15 The Magnus moment (Lieske, 1966). 

D. ESTIMATE OF REQUIRED CHANGE OF MAGNUS COEFFICIENT 

The desired change in the Magnus moment coefficient can be 
estimated by varying the values of the coefficient and plotting the 
results as a function of S 2 , S 3 . This is done in Figures 16 and 
17 below. The first Figure shows how the terms have a polarizing 
effect since they contain opposite signs in S 2 , S 3 . This is 
desirable since any change in the value of the coefficient may 
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cause an instability to result from either term. The second plot 
is the product of S 2 > S 3 demonstrating the narrow region where 
stability is satisfied (the region where the product is positive). 
For typical values of the coefficient, a magnitude change by a 
factor of five may be sufficient to cause instability. 


MAGNUS MOMENT COEFFICIENT vs.STABILITY TERMS S2 and S3 



-0.5 0 0.5 


Magnus Moment Coefficient 


Figure 16 Stability factors vs Magnus moment coefficient value 






























S2 x S3 


THE S2 x S3 PRODUCT 
vs. 

MAGNUS MOMENT COEFFICIENT 



Figure 17 The stability factor product produces a region of 
stability for certain Magnus moment coefficients. 

E. ESTIMATING YAW CHANGE VIA CHANGES IN THE MAGNUS MOMENT 

The next set of calculations conducted to provide a theory for 
instability are based on a point mass trajectory model developed by 
the Ballistic Research Laboratory. The equations derived for this 
model provide a simplified means to calculate the yaw of repose. 
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a e , a dimensionless expression to quantify the effect of changing 
yaw angle along the trajectory. The goal of the BRL was to develop 
a model that would incorporate the effects of yaw without the 
computing time of a complex rigid body simulation. The estimate 
was an intermediate step in achieving this goal. For this 
discussion, 3 e is the desired term for examining a magnitude effect 
of varying the Magnus moment (Lieske, 1966). 

If the results of this calculation demonstrate a dependence of 
the yaw of repose on the Magnus moment, then the stability argument 
previously suggested is validated. The BRL found that the 3 e term 
was successful in matching empirical data when applied to equations 
of motion that included yaw drag and lift. They did not apply this 
term independently to varying conditions of the surface roughness 
of the projectile. 

In the development of the yaw of repose estimate, the BRL used 
some notation that is not consistent with previous chapters. To 
avoid confusion and define some new terms, a list of these symbols 
is provided below (although Lieske et al. , accounted for the 
Coriolis force, it has been discounted by many other studies as 
negligible and is ignored in this magnitude analysis as well): 

B Transverse moment of inertia 

g Q Acceleration due to gravity at ground level 

T Unit vector in the direction of v 

K n Drag force coefficient 

o 

K n Yaw drag force coefficient 

u a 

The same coordinate system previously introduced, a ground- 
fixed right handed Cartesian system, will be used. The body is 
still considered a solid of revolution and the unit vector x is 
along the axis of symmetry pointing forward through the nose. The 
magnitude of the angular velocity, N, is parallel to 51 and is 
positive for a right-handed screw motion along x. 

The following are the simultaneous differential equations of 
motion for a spin stabilized projectile. First, for the center of 
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mass: 


Pd 2 

m 


(K + & 2 K d ) 

^ O OL 



V* (x x v) 


pd\. -r , Pd 3 ^ 

-—K vx +- 

m s m 


(x*v) +g. 


(3.7) 


Recall that the total angular momentum H can be expressed by: 

H = ANx+B (x*5) . (3.8) 

Therefore, the sum of the applied moments gives the vector 
rate of change of the angular momentum: 


fi=ANx+ANX+B{x*X) =p d 3 K M v(v*x) -pd^K^ix^X+pd^K^lx^ (x»v] -pd i K A Nvx. 


(3.9) 


Defining the yaw of repose from the rigid body system of 
differential equations requires these assumptions: 

1) The projectile can be represented sufficiently as a 
body of revolution. 

2) The projectile is dynamically stable. 

3) The initial yaw is assumed to be small and has a 
negligible effect on the trajectory. 

In equations 3.7 and 3.9 the magnitude Jx* 171 is present in the 
lift and Magnus terms. If we let I=—, then the magnitude 
Ixxj|=sin5. Lieske, et al., then define the yaw of repose: 

3 e sfx (xxi) =x- C os5J. (3.10) 

This vector represents the vector yaw directed from f toward J? and 
is in the plane determined by I, x. The effect on the trajectory 
by a e will be small with the assumed conditions. Therefore, (E e 
is assumed negligible yielding the following approximations: 
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5 = cos 8i, 5 = cos8j. (3.11) 

Then, dividing the total angular momentum into parallel and 
perpendicular components results in: 

AN=-pd i K f Nv, ANX+B{x*X) =pd 3 K M v(v*x) -pd 4 K H v(x*.*) +pd i K T N[x* (x*v) ] . 


Then substituting x and its derivatives in 3.7 and 3.12 with 
equations 3.10 and 3.11 will provide a solution: 


u= 


-pd 2 {K D +b 2 K D ) v 2 I pd 2 K 1 y 2 (i e pd 3 K s vcosbl pd 3 K p Nv 


m 


m 


m 


m 


5/1) +9- 


(3.13) 


ANcosfi J +Bco s 6 +Bcos 2 8 (f»?) =pd 3 JC M v 2 (f«a/ +pd 4 JC r Wv[cos8 (a e +cos6l) -I] -pd 4 K tf v r cos6 [ (a^+cosCJ) *1] . 


Cross multiplying both sides of equations 3.13 with I and utilizing 
the terms, I, I the yaw of repose becomes: 


[-AKjNcosft ( V*v) +md 2 K J JHcos6 [?* {u-§+pd i K 3 cos6^ r /m) ] -K L [v- (v'T) T] [2Bcos 3 5 {vT) -pd*K H cos 2 6v 3 ] ] 
‘ [pd 3 K I iC^v 4 + pd 5 J<' y iiC 3 jW 3 cos5v 3 +iC L Bcos5 [ (v*f) 3 - (v-v) ] ] 


(3.14) 


With the v term dominating the denominator, the last term can be 

2 

discarded, leaving a ratio of . The middle term in the 

v 4 

denominator after Kt can be considered approximately equal to g. 
With some known behavior of artillery shells this further reduces 
to a final form: 

-AKjJsKvxti) +md KyN[v* (3-g) ] (3.15) 

S =- . 

0 O yl p- 

Pd K L K M v + P d KpKjN v 


The next step is to use this equation to estimate the effects 
of varying the K terms. v*^ is estimated using the ideal 
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2 . 2 * 

trajectory, this is done by approximating a e =sm 0 
can be written: 


Then, 3.15 


c£ e = (a b -a a ) (v*u) -ab(v*§) , 


(3.16: 


such that, S e « (a a -2a^) (v Q *§) 


(3.17! 


Now the useful form desired for sin 5 can be written: 


2c 22,2 2 . 

sin 8 = v Q g (a -+411^-4 a a a^) 


13.18) 


For the mathematical analysis it is necessary to define the 


terms oc a , . 


a * = —3 


ak l n 


[3.19) 


p d 3 + p d^KpK^v 2 


mKjN 

^ jb . o o 

pdJC L %v +pd iCpK-yN v 


It is worth noting that again, both terms contain coefficients 
for the Magnus moment. Using numerical values, and changing the 
Magnus moment coefficient from 0.3 to 1.5, the solutions for the 
equations 3.19 are: 


a a « 9x10 


-9 sec' 


~ 4x10 


-8 sec' 


Thus the result for initial values of yaw will be: 


oJ a ~ ixio 5 * sin 2 5, 8 ~ /TOxio 2 . 

Similarly, if the Magnus moment is increased by a factor of five, 
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then the results are: 


a a ~ 9»10 9 , a b ~ 7xio 7 , Sg * lxio 3 , b ~ /TOxio 1 . 

This is an order of magnitude change! The theoretical evidence to 
support instability with a five-fold magnitude increase in the 
Magnus moment is established. These results although not linked 
directly to the stability criteria support the argument that 
instability can be achieved with a change in the Magnus moment. 
Based on the stability equations, we expected instability to arise 
by changing the Magnus moment. Here, we have changed the Magnus 
moment and verified that the yaw angle will change considerably. 
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F. 


WILL THE MAGNUS EFFECT CHANGE SUFFICIENTLY WITH 
INCREASING THE SURFACE ROUGHNESS 


A drawing of the projectile as viewed from the nose tip coming 
out of the page will be helpful for estimating the magnitude of the 
Magnus effect. 



The cross flow, v ± , is approximately 50 m/s for a small yaw of 
about five degrees. This is nearly equal to the velocity of the 
boundary layer close to the projectile surface (laminar flow) which 
is calculated to be 60 m/s. Applying the Bernoulli equation to the 
unroughened laminar flow projectile surface leads to the following 
analysis: 



= constant = P Q + 



(3.20) 


where v-^v^ giving 
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4] 




Adding a surface roughness, as in Figure 19, to the entire 
body would cause the resultant velocity vector on the lower side to 
decrease as fluid momentum from outside the boundary layer is drawn 
into the boundary layer by turbulence. The velocity along the top 
of the projectile will not change appreciably since the component 
of cross wind velocity and the velocity of the boundary layer 
(caused by the angular velocity of the spinning projectile) are 
nearly equal. The result is a small change in the Magnus force 
that will be estimated now. 

For the projectile with a symmetric band of roughened surface, 
the view from the nose shows how the turbulence draws the 
surrounding flow into the boundary layer in the roughened region. 
The turbulence transports fluid momentum across the boundary layer 
and tends to bring the flow velocity closer to the surface. Notice 
the turbulence does not change the profile of the velocity along 
the top surface of the projectile. 

The turbulent boundary layer is turbulent at a Reynolds 
number, 

„ . p v- 


in a band around the projectile, using the full flow velocity 
including the longitudinal component because it is the full 
velocity at the boundary layer edge that determines the Reynolds 
number. 
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unchanged profile on top 



Then comparing the equations in 3.20 results in the following: 


P 2 + |pv 2 = P 0 + |p v? 


(3.21) 


approximating v i ss "x v j_' then 


1.21-2. _ 132 

p 2 * P o* 2 P <’'x- 4 V i> * P 0 + 2 P 4 V J- 


such that now, 

A p turb “ "g P v x ' "* ^ p laminar~^ p tuzb ~ ’ 

This is a change from the laminar case by a reduction of about 25% 
change in pressure and hence a similar change in force. But, this 
is only applied to the band of roughness which for a laser induced 
case, may only be one tenth of the length of the projectile. As a 
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moment arm, this would represent 20% of the moment arm length. 
Thus, applying a 25% decrease in Magnus force to only 20% of the 
surface would result in about a 5% change in the Magnus moment. 
This is not a significant enough effect to achieve the instability 
we had hoped for via the Magnus moment. 

If the roughness were applied to the entire surface, the net 
result would only be a change in the Magnus force by a total 
magnitude of about 25%. This is still far short of the five-fold 
increase which is likely to cause an instability. There may be 
additional non-linear factors that contribute to the Magnus effect 
such as those described by the vortices that produce an additional 
sideforce (Platou, 1956). The analysis applying this additional 
component may still provide validity to the prospect of exploiting 
the Magnus moment as a means to achieve instability, but is not 
undertaken here. 

The case of an asymmetric roughness does not present a uniform 
force that could be applied to the Magnus effect analysis. In fact, 
the asymmetry may only cause a wobble of the same proportions 
derived above. The asymmetry may also appear symmetrical to the 
Magnus effect if the spinning boundary layer creates a uniform 
turbulence at high spin rates. However, it may create a moment 
arm of drag which will be discussed later. 
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IV. AERODYNAMIC DRAG AS A MODEL 


The drag will be examined as the friction factor changes by a 
factor of five; this value was not chosen arbitrarily, however the 
analysis of the Magnus effects proved that a similar change in drag 
would not result in an instability. As such, it is hoped that a 
similar magnitude change in the drag will result in an additional 
drag force sufficient to reduce the range of the projectile. Drag 
will be examined separately to determine if the effects of drag 
alone can prevent an incoming ballistic projectile from reaching 
its target. 

The model for the analysis will be a cylinder to approximate 
a projectile. The flow through a cylinder has been, in particular, 
very well studied. A pipe is a cylinder with internal flow that 
approximates the flow behavior on the external surface as well. It 
has been well established that increasing surface roughness in a 
turbulent flow will have an effect on friction resistance. In this 
case, for the projectile, the effect is on the aerodynamic drag 
across the whole body. For the model, then, we want to increase 
friction resistance by a substantial factor of three to five. 

The friction resistance is due to the surface effect of 
changing the boundary layer. By increasing surface roughness, the 
sublayer will break up and become increasingly turbulent. 

For the 155mm ballistic projectile a roughness value, e, of 
.046mm (White, 1979) can be divided into the length (700 mm is used 
as an estimate to account for the geometry) to enter Figure 21 for 

friction with smooth and rough walls. With a value — of 15000 

7 e 

and using a typical Reynolds number, 10 for the projectile, the 
flow can be determined as past the transition point, but not fully 
turbulent. Therefore, the condition of turbulence is met and an 
increase in surface roughness should result in dramatic changes in 
the drag. In the next section, the calculations will be made that 
show the actual values of the Reynolds number and how the roughness 
must vary to provide the necessary results. 
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Figure 21 Reynolds number vs drag coefficient 
(White, 1979). 


ft.. AERODYNAMIC DRAG ESTIMATES 


The Reynolds number for our projectile must first be 


estimated: 


Bk « 2 . 6 * 10 7 


(4.1) 


■L V 


Then, entering Figure 21 with that result, and moving acr< 




the chart with the value of 


= 15000 to determine the drag 


coefficient, gives C D = .004. An increase of more than three can 
be seen on this Figure using C D £ .013 resulting in an of 
approximately 200. This results in a required roughness of 3.5 mm, 


€ = vrvv = 3.5mm (4.2) 

200 

Similarly, if the drag coefficient is used as the measure to 
determine a change in roughness required, it is possible to 
calculate the coefficient by (White, 1979) : 


k d = 


■ 031 
R 1/7 


8700 


(4.3) 


Using this result and multiplying it by five will again result in 
a change of roughness L/e on the order of several hundred and a 
resulting roughness of several millimeters. This roughness can be 
visualized as dimples covering the surface of the projectile that 
are several millimeters deep and spaced so that no smooth surface 
is larger than the diameter of those dimples. The number of 
dimples required can be calculated using the area of the projectile 
and a dimple size of 5 square centimeters to be approximately 
162,000. This means that the laser must pulse that many times on 
the incoming projectile. With a short pulse length on the order of 
nanoseconds, this can be done in milliseconds. A longer pulse 
length but no longer than microseconds will still allow sufficient 
roughening in under a second. With the coning and spinning motion 
of the projectile, a pulse rate somewhere in the middle will allow 
for the most dispersion of pulses over the entire body. 

Similarly, a continuous pulse laser may melt the surface 
sufficiently in less than 1/2 second to cause the same effect. The 
roughness does not need to be caused by dimples. It may be equally 
effective to produce ridges of the same magnitude of solidified 
molten flow produced by the laser energy for a roughening effect. 

From the criteria for vaporization equation, an approximation 
of the energy required can be made (Schriempf, 1974) : 
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F 0 X p ~ P& x (C s {Tm Tq) +l M' F o x p ~ (56 00 -), (4.4) 

cm 

where Ax is the hole size or roughness desired. Substituting in 
the Ax = 0.35 cm term results in the following requirement 
(details are located in Appendix B): 



The symmetry of the roughness does not significantly affect 
the net force acting on the projectile. Because the projectile is 
spinning and the roughness will be asymmetric over the entire body, 
the turbulent boundary layer will spin with the projectile as the 
cross wind interacts with the turbulence. In fact, an asymmetric 
disturbance of several thousand holes around the circumference of 
the spinning projectile will produce the largest increase in the 
drag effect since the boundary layer and hence surface friction 
will increase around the entire surface. 

B. DRAG VS. RANGE ESTIMATE 

The question remains, "how does this drag reduce the range of 
the projectile?" To complete this estimate, the equations of 
motion incorporating the drag will be used. The vertical drag will 
be assumed negligible compared to the axial drag for ease of 
estimation. This is probably a reasonable assumption anyway since 
the projectile will be engaged close to its terminal point and the 
altitude will be small compared to the range. Additionally, the 
yaw as the projectile reaches its target may be at its maxima 
providing additional lift to counter any vertical drag. This 
assumption also allows a starting calculation for the time of 
projectile flight which will be needed in the more detailed 
equation of motion. 
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v y t- —gt . 


(4.5) 


Using a firing angle of 45 degrees for the projectile to achieve 
maximum range, the horizontal and vertical components of the 
velocity are equal initially ( V v =V Y ) . The projectile is 
calculated to have 79.5 seconds as a time of flight. 

The equation of motion-used to determine range information is: 


x = - - p (Y) - j-r-xKp = -G(v) x. 

m u 5 (Yj v 


The functions that, in general, are dependent upon altitude will 
also be considered constant for the estimate made here. Applying 
the following steps to equation 4.6 results in a final form useful 
for determining the range shortfall with an increase in drag: 


Finally, 


so that 


x = u, u = -u G(v) , In U \t) = -g(v) t, 

u(o) 


u(t) = e -u(o)G(v)t_ x(t) = j- e -G(v) t dttXo _ 


-G(v) t. 


x( t) 


v _ V 

x o -G{v) t „ u r, X ° 
- e +X., where X=- - 

G(v) ° ° G(v) 


o . ^ -G(v) t. 


x(t) = —If (1 - e 
G(v) 


Applying the final form of 4.7 to the projectile's flight in 
two steps will yield the desired results. First, by choosing a 
point 2000 meters from the CBWS the projectile's time of flight and 
range up to that point can be calculated. Then, using that point 
as the start of the engagement with a the new drag coefficient, the 
range of the engagement leg of flight can be compared to the range 
if the drag were not changed. The result, for our increase in drag 
is a 20% reduction in range with the engagement starting at 2000 m. 

From the launch of the projectile until 2000 meters away from 
the target: 
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For the remaining 2000 meters of travel new velocities in both 
the vertical and horizontal components must be applied in addition 
to an increased drag coefficient (increased from .004 to .02). The 
velocity will be calculated from: 

x(t) = v x e t (4.8) 


where, G= .002283 for a drag factor of .004. This will give an 
additional projectile range during the increased drag portion 
(G= .008121) of its flight path of only 1619 meters. The 
projectile will fall short of the target by nearly 400 meters or 20 
percent of the engagement distance. A twenty percent reduction in 
range with an easily obtainable laser requirement is an encouraging 
result! 

C. YAW TORQUE CREATED BY A CHANGE IN ROUGHNESS 

In the case where only one pulse or a series of pulses in 
close proximity roughen the surface of a projectile in a small area 
asymmetrically, there will be a torque induced that will produce an 
increase in yaw with the same frequency as the spinning projectile. 
Depending on the existing yaw angle and the cone rate, this torque 
may cause tricyclic motion of the projectile. 

If the asymmetry is out of phase with the cone rate, it will 
cause a coning motion within the larger coning motion, thus a 
tricyclic motion as in Figure 22. This will almost always be the 
case. There will be points in the tricyclic motion where the yaws 
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are complimentary and result in a larger yaw and points where the 
yaw motions cancel. The effects of this complex motion are 
difficult to describe and would vary in almost every case. 

With this in mind, the projectile will be assumed to have zero 
cone for this calculation initially. This will enable a 
calculation of the contribution of the roughness independent of the 
other forces. 

Dividing the drag force into components, the component along 
the axis will cause a yaw to develop. The torque created by this 
force is solved by the following equations: 

T = a*F = 131 lFlsin<l> (4.9) 


51 






= I SI IFl sin( 180 -P) 

* 151lFlsin(180-Y) 
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a 

« r*F, where r=7 cm « R . 



Figure 23 Asymmetric drag. 


An equivalent expression is then used to solve the equation in 
terms of flight time remaining to the target. At. The change in 
angular momentum is related to a change of yaw angle that is chosen 
large enough to increase drag sufficiently to miss the target. The 
result is expressed below as a sequence of steps assuming there are 
three seconds of flight time remaining: 


F-r 


Al . L s 0)Att 
At At 


where L s = 



(4.10) 


This results in a required force of approximately 300 Nt. But, 
this is not a useful number unless we can attain it with our 
estimated roughening that increased drag by a factor of five. For 

this, 

F = —pU 2 C D -A, SO AF = pu 2 (A C d )A « 35 Nt. (4.11) 

2 u * 


This result is an order of magnitude less than the 300 Nt required. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The stated goal of attaining a system that will perform the 
requirements set forth in the Defense Planning Guide is within the 
realm of physical possibilities. The use of directed energy may be 
a viable solution to the problem of defeating an incoming spinning 
ballistic projectile. 

The following potential mechanisms were investigated with 
available theoretical models and quantitative estimates were made 
with the following results: 

(1) Impulse transfer of momentum of 10% of the projectile 
momentum by laser ablation: The allowable laser energy fluences 
before laser absorption waves will decouple the laser from the 
target are orders of magnitude too small to achieve the required 
lateral momentum change. 

(2) Impulse angular momentum change by laser impulse to 
increase the yaw angle significantly to induce large drag: Because 
of the limitation of laser pulse energy by the absorption wave 
threshold, it would require multiple sub-nano second multi-gigawatt 
pulses to achieve significant angular momentum changes. 

(3) The stability of the projectile and its dependence on 
various aerodynamic coefficients was investigated using a well 
known stability model. The Magnus force and Magnus moment were 
identified as the most likely candidates for affecting 
destabilization by a five fold increase of the Magnus force. It 
was shown that such an increase would indeed lead to a significant 
change in the yaw. However, a simple analysis of the effect of 
surface roughening by a factor of 3-5 by laser melting or dimpling 
leads only to a 25% change in the Magnus force because the 
turbulence tends to wash out the pressure differential between the 
top and bottom of the projectile. 
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(4) The effect of increasing aerodynamic drag by laser 
roughening of the projectile surface was investigated. It was 
found that a 3-5 fold increase of surface roughness (drag) would 
require fluences of 2-5 Kilojoules per square centimeter and that 
the effect of the increased drag could lead to a 20% decrease in 
the projectile range over the last 2000 meters. 

(5) Creation of a yaw increasing torque by a single 
asymmetric high roughness laser melting spot on the projectile was 
investigated and it was found that with a fivefold increase of 
roughness on a 10 square centimeter spot, the torque is too small 
to cause significant yaw increase and subsequent drag increase and 
range reduction in the time available for engagement. 

The aerodynamic drag was certain to increase regardless of the 
Magnus effects when the surface of the projectile was roughened. 
The result although not dramatic emphasizes the viability of using 
a directed energy CBWS as a means to achieve a "kill" of an 
incoming enemy projectile. If the target of the projectile were a 
ship, then the range reduction would certainly be sufficient. For 
land targets, the decrease in range may pose other problems in a 
congested or civil warfighting environment. 

The system requirements may be of a scope that is entirely 
different in nature. It is necessary to develop a system that is 
capable of tracking a projectile for ranges of several hundred 
meters during an engagement. This must be done with an accuracy of 
centimeters. The system must also be capable of determining the 
exact location of the projectile relative to the projectile's 
intended target. It would not serve its purpose to engage a 
ballistic projectile that was certain to overshoot its aim point if 
doing so would cause it to fall short and on target. These are 
just some of the system problems that will require addressing 
before a suitable CBWS could be built. 

The energy requirements for the laser to induce sufficient 
drag are well within the bounds of current technology. In that 
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light, it is advisable to find an experimental means of validating 
these findings. There may be some interest in determining exactly 
how deep the laser effects will penetrate or what size melt ridges 
will develop on the ballistic projectile surface when atmospheric 
disturbances have reduced the effectiveness of the laser. A sea 
environment in particular, may have a deleterious effect on the 
directed energy against a relatively low flying target such as a 
terminal ballistic projectile. 

Other recommendations for research include a more detailed 
look into the drag behavior on a projectile. The model of a flat 
cylinder does not take into account such things as the point where 
the shock wave detaches from the body. There may be other unique 
aerodynamic behaviors with various projectile body shapes that 
would require significantly more or less roughening to achieve the 
desired results. 

Finally, one last attempt to exploit the Magnus effects is 
worth the effort. The nonlinear effects of the vortices shedding 
may be sufficient to produce the desired instability. It may be 
more desirable to pursue the use of a kinetic energy vehicle to 
cause a change in the center of mass location such that the Magnus 
moment causes instability. This would likely be a solution that 
would provide a much larger safety margin of destruction from the 
projectile's intended target. 
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APPENDIX A 


Aerodynamic Forces Acting on the Projectile Preliminaries 



CP s center of pressure 
R = resulting aerodynamic force 



Figure 24 A single resultant 
force on the projectile. 




For this introduction into the stability a preliminary model 
will be used to define some of the most basic motions and forces 
acting on the projectile. The projectile in the Figure above is 
rotationally symmetric with velocity, u , with respect to the air 
and an angular velocity, (0 , equal to zero. The terms of this 
Figure are defined below: 

M = Overturning moment = component of T along the Z-axis 

D = Axial Drag = component of ft along the X-axis 

L = Lift Force = component of 5 along the Y-axis 

M, D, and L depend on the yaw, 6, by: 

L=X sin6 X, m are independent of 6 

M=msin8cos8 

D=pu 2 d 2 K D , L=pu 2 d 2 X L sin8 , M=pu 2 d 3 X M sin8cos8 
X=pu 2 d 2 K L , m=pu 2 d 3 K M 

D, m, and x are all functions of the projectile shape, 
diameter, velocity, sound speed and density. 




a 


Figure 25 The forces and 
coordinates. 
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Since it is desirable to have the projectile travel with some 
yaw, it is necessary to impart a stabilizing spin. This will lead 
to a complete description of the forces acting on the projectile by 
reducing the aerodynamic forces to a finite number of coefficients. 
These have been verified experimentally within a satisfactory 
degree of accuracy for describing the motion of the projectile with 
a small yaw. 


?((3,u)= net aerodynamic forces acting at the center of mass, 
3 ((S, u) = net torque due to aerodynamic moments. 

,% 2 /^3 = unit vectors fixed to the projectile. 


F=F i ^ i 
G=G i J? i 
U=U £ 

CO =(D 


CM. 


Then solving for F£ and G^ is done with a Taylor expansion 
about the unyawed position. The coefficients of the expansion are 
independent of u and (0 to the first order. Adding s, the angle 
of rotational symmetry between zero and pi, the complex numbers, 
£=u 2 + iu 3 ; T) = o > 2 +iw 3 where u 2 =— (£+0 ; u 3 =—i($~ 5 )/ become 
necessary. Linear functions of u 2 , u 3 , <0 2 , <<>3 can be written as 
linear functions of £,£,T|,T). The Taylor expansion then becomes: 


Ff = a i +b^ + b 2 T| + + jB 2 t[ 

^ = F 2 +iF 3 = a 2 +0-^ + c 2 ii + d-L^ + d 2 rf 


(A.l) 


The coefficients except a are in general complex and the special 
form of F is necessary because it must be real for all values of 
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£,r| . Rotating the vectors ,& 2 >^3 through the angle s about the 
axis of symmetry changes the axis system to >&2 '^3 an< ^ t ^ rie 

following transformations take place, Z ~ Z e 1S , T) - r| e 1S ,Z - Z e 1S 
rf - r| e , F ± - F 1 , and 3F e . Because of the symmetry, such 
a rotation cannot change the coefficients in equations A.l. The 
resulting equations are: 


a l + 


b ± Z e is + 


b 2 r) 


e~ ls + 


b ± Ze 1S + b 2 ne 1S 


(A.2) 


e ls = a 2 +c 1 £ e is + c 2 x] e + d^e + d 2 rie 


Comparing equations A.l and A.2, the terms on the right from before 
and after the transformation can be set equal in F 1 for all 5, T| . 
Choose r|=0, Z= ^ -* 2b 1 b 1 =2b 1 ^coss, coss^l - b-^0, similarly b 2 = 0. 
Using the same equations to compare &, the d terms can be examined 
the same way and they result in coefficient values of zero as well. 
Similar arguments can be made for the components of G. The result 
is that if a projectile is symmetric about its axis of rotation, 
its force and momentum system can be written to a first order 
approximation: 

F 1 =a 1 &~=F 2 +iF 3 = C 1 Z + C 2 r \ 


(A.3) 


G 1 -e 1 &=G 2 +iG 2 = C 3 £+C 4 n 


where a-^,e-^,C-^ ...c 4 are functions of density, shape, velocity, 
angular velocity, and wind speed and are independent of the 
orientation of the axes ,St 2 ,£3 with respect to the body. One 
can further show that if the projectile has a plane of mirror 
symmetry, then in A. 3 the terms F ± , Stc-L , 9tc 4 , %c 2 , Sc 3 are even 
functions of and ^’ 1 >^ c 2 , ^ c 3 '^ c l , ^ c 4 are odd functions. 
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The resulting dimensionless form of the force and moment 
coefficients are: 

C -l = -pd 2 u-L%+ ipd 3 co 1 iC F 

4 . 3 

C*2 = p C? G) 1 ^JKF ^ P ^ 

=:— pC? (D-^.FCiji'-.Lpc/ (A- • 4 ) 

4 .5 

“ p C? "** P ^ ^l^XT 

*1 = -pd 2 UiF DA 
Gi = -pd 4 u-LF A . 

The K 's are all even functions of &)-, and of —, d— . The signs 

x a u 

are chosen so that for a typical projectile all coefficients would 
be positive (for most shells, K T is negative) . 

The forces and moments are then summarized as: 

2 

pd u-> KJc, The Normal Force cross force due to cross 

velocity. 

3 

ipd The Magnus Force Magnus cross force due to 

cross velocity. 


pd 

ipd 3 u 1 F 5 ti 
p d CO K j,£ 
ipd 3 u 1 F M ^ 
pd u-j 


The Magnus Moment 
Over Turning Moment 


Magnus cross force due to 
cross spin. 

Cross force due to cross 
spin. 

Magnus cross torque due to 
cross velocity, 
cross torque due to cross 
velocity. 

cross torque due to cross 
spin 


Damping Moment 
or Yawing Moment 


5 

i P d (x) 

_ ,2 2 

pd u x iC DA 

4 

pd u 1 co 1 K' a 


Magnus cross torque due to 
cross spin. 

Axial Drag 

Spin Decelerating Moment 


The force R has been resolved into forces along and 

perpendicular to the projectile instead of u. Taking into account 

2 

the yaw, 8, the normal force becomes pd (ucosfi) X^-(usin8) , and the 
axial drag becomes pd (ucos5) K DA . Computing these with 
components along and perpendicular to the velocity vector gives: 


pd z u z K L sin6=pd z u z K N sin5cos 8-pd u Jf DA sin8cos^6 


pd 2 u 2 K D =pd 2 u 2 K N sin 2 5cosb + pd z u z K DA cos^6 . 


(a. 5 ; 


A comparison then gives: 


2 2 
k l =k n° os ^ "“TC^^COS ^ 


(A.e: 


2 3 

k D =k N sin ^cos6+k da cos 8 


which to a first order approximation (8<1) is: 

K L* K N~ K DA an< ^ K D~ K DA • (A. 7) 

While the moment term, K M , remains identical to the one previously 
introduced. 

The general equation of motion for a symmetrical projectile 
has four assumptions: 
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1) A fixed earth three dimensional coordinate system ( X, Y, Z) 

neglects the Coriolis force. 

2) Values such as density that change with altitude are 

assumed constant. 

3) There is no wind velocity. 

4) The projectile has an angle, s, between zero and pi of 

rotational symmetry and a plane of mirror symmetry. 

The velocity vector, angular velocity vector and projectile 
fixed coordinate system remain the same for the notation of this 
section. Equations A.3 represent the Aerodynamic force and the 
Aerodynamic torque, respectively. A unit vector parallel to the 
vertical axis in the earth fixed system is defined to be 
Y-]_ Y=Y' 2 + ^ 1’3 and the acceleration due to gravity transforms 

~gY 1 - -gT. The remaining terms to introduce are: the moment of 
inertia about the axis of the projectile. A; the moment of inertia 
about any axis perpendicular to the projectile axis through the 
center of mass, B; and the mass of the projectile, m. 

The general equation of motion is: 

^(m,u) =F-mgf. (A . 8) 

The angular momentum equation in the projectile coordinate 
system is given by: 

2! =AC0 1 3? 1 + B(jb2&2 + SM 3^3 = (A~B) 0) 1 ^ 1 +B(5 

(A.9) 

and —— = G , —— { (A-B) G).. & +BC&} = G. 

dt dt 11 

The time derivative of vector P is expressed in the projectile 
system where P transforms, P- {p^; P = P 2 + iP 3 ). Using index 
notation, P=P^^ and P=Pj£j+PjX. The vectors £j are fixed in 
the projectile. 


P = Pj£j+PjG>x£j 
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Xj = (x>x£j 


(A.10) 


P'&k ^jk + Pj ‘&k ~ p k p j ^ ^k x ^j^ 


In components, the first one becomes: 


P -£i =P 1 + [ (P 2 +iP 3 ) («2' <0 3 ) -(P 2 -iP 3 ) ( w 2 + i(0 3 ) ^ T =P 1 + ( p ^- pT l) 


(A.11) 

The second and third added together with the third multiplied by i 
becomes: 

P-£ 2 + iP-£ 3 = P 2 +iP 3 -iP l (<«> 2 + iG> 3 ) +ito l ( p 2 +iP 3^ = P-iPi^l +i(0 1 P ( A .i2) 


So, the equation of motion in projectile based coordinates are: 
m[u 1 -(^ri-fTi) -|] =F 1 -mg 1 


i[$-iu 1 T|+i<D 1 5] =&'-mgT = C 1 i;+C 2 'n -mgT. 


(a. 13 ; 


And, the equations for the angular momentum are: 


( A-B) co 1 ^ 1 +B(3 =G. 


(A.14) 


The above lema on the time derivative gives 


— [(O-L^] ={ <b 1 +(0-Ti-0-Tl) -2 = <J>i 
{ 0 -iw-lTI + 0 = -io> 
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d d. 

— [o5] ={ <!) 1 +(Ti*ti-ri-'n) - = 6 1 

{ f| +ior| = f| . 

With this we have the two equations: 

( A~B ) + = 

(A.15) 

(A-B) T|) + B-i ) = r = C 3 Z + C 4 T1 


The differential equation for a change in the force of gravity in 
the projectile frame is given by assuming Y = constant, Y = 0 in the 
lab frame and applying the lema: 

4z [?] = { Y, + (Tii-TTi)4 (A. 16) 

at x 2 

{ T - iY 1 1| + icoT. 

There is however an altitude dependence on the speed of sound 
and the density. The altitude change of the projectile must be 
tracked in order to insert the correct p(l r ), a (Y) values. If u has 
a component ft? in altitude direction then its altitude change is 
given by: 

y = ft? = u-l^ + u 2 y 2 + u 3 y 3 = 12^ + ~ (S? + IT) . (A. 17) 

The normal equation of motion is then derived. Again, several 
simplifications are introduced to get an approximate trajectory 
before making corrections: 

1) The projectile axis is tangent to the trajectory. 

2) The only forces acting on the projectile are drag and 

gravity. 

3) The variables affected by altitude are considered constant. 

4) There is no wind. 

The projectile will remain in the vertical plane of the initial 
velocity vector as indicated by Figure 26. 
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Figure 26 The drag opposes the motion. 


D = pd 2 U 2 K D , u 2 = X 2 + Y 2 , D - pd U Fg = 

The equations of motion are then: 

mi = -pD uK d u - mgf (A. 18) 

with components: 

irdt = -pd Z uK d X, mt =• -pd 2 uK D t - mg (A. 19) 

where K D = , u 2 - X 2 *? 2 , P=P (« , a-a(r> . Most exterior 

ballistics involves the precise integration of equations A.19. 

The standard form of a normal solution is given by. 

X = ucos0 , Y = usin0 

and X = ucos0 - usin0 0 , Y = usin0 + ucos0 0 , so that 

imicosO - musin0 0 = -pd 2 K D u 2 cosQ (A.20) 

2 2 . 

and musinO + mucos0 0 = -pd XpU sin0 - mg. 

Multiplying the third and fourth equations above by cos and sin 
respectively, then adding the result gives: 





(A.21) 


2 2 

mu = -pd K D u - mgsin0. 

Similarly, multiplying by sin and cos respectively then subtracting 
gives: 


0 = - —COS0 
u 


(A.22) 


To the above equations A.21 and A.22 the following terms are 
added: p = p(y) , a=a(Y) , K D =K D (~) . 

Next, the angular motion of the projectile will be discussed 
starting with the simplification of the equation of motion, 
equations A.13 and the following: 

Aob ^ = G-^ 

Bf 1 + (B-A) iO) 1 T| _= g 

Y ± = (Tri - ft]) ± (A. 23) 

T = iY 1 n - ico-i^T 

* | (IT * IT) 

where = C^ + C 2 T), % = C 35 + C^r\ , and the values . . .C 4 , F-^, G- 
in terms of the aerodynamic coefficients are given. 

•, Changing the notation somewhat to include a density factor, 

p-, is necessary. The density factor is applied to •, the 

m . . d J 

aerodynamic coefficient to result in J terms, J-. = p- K n . 

m cP 

Further, we will define the spin per caliber travel as: v = < 0 -,-. 

U 1 

The coefficients then become: 

U 1 

- ( Gfj + ivJp)m-^j- 

C 2 = (v J ^ 1 G^) mu^ 

0 2 = (“"V — iJ^jmu^ (A. 24) 

~ ( — Gjj ^ inu^ d 

G-^ = ~J m • 

Then use approximation by means of the normal equation of 
motion. Let the center of mass motion be governed by the normal 
equation of motion: 
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■pd 2 u 2 -^ - gsin0 = -u 2 -^ - gsin0 


-gc os- 


(A.25] 


X = ucos0, Y = usin0, 


The solution is assumed to be known and satisfies 


U = -Jp-T ~ gsin0, 


(A.26) 


Next, we will conduct a change of independent variables 
starting in the yaw equations. First replace: 


r U 1 r U 

p ’f ~d dt -f { d )dt - 


(A.27) 


This is then the arc length along the trajectory in calibers. The 
equations determining the yaw motion will become essentially 
independent of the size of projectile. A large projectile has the 
same period of yaw measured in calibers of travel as a small model. 


(f» - =- Q' — 

dP d dP d u d 


(A.28) 


The trajectory normal equations are then: 

0 0 

u' = -J D u ~ gdsin-y., 0' = -gdcos— 

U jjZ 


(A.29) 


The equation for the derivative of the angular velocity becomes 


Ati) ~ GCO 


dG -i J a Tj n v 

— = -dvm— 4- with o, - -^v, co' = v'- + -U' 
AU A U 1 d d d 
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u V U V 

Therefore, v — + —-C7' = - dvmJ A — , v ' = -—[/' 
rl n /i 2 i rr 


2 ^A 

d vm-^- resulting 


v' = -—^md 2 + \J D + vgdd 2 sin 0 . 


(A.30) 


To change dependent variables the variables, 5, r|, , T are 

changed as follows: 


A- I. 0 


L f vdp 


J. 

if vdp 


y = gTe 


d 0 


gi = gr 1 


/vdp 


— is the vector yaw which has the magnitude of the sine of the 
angle of yaw. The factor 

P t 

eiifvdp] =e[iJco 1 dt] 

0 0 

changes the coordinate system from ^ / £3 which rotates with 

the projectile to one that does not rotate about the j^-axis. 

The quantity, A, , is the complex yaw measured from the axis of 
the projectile to the trajectory on a coordinate system with one 
axis along the projectile axis and which is not spinning about the 
axis. Similarly, |i, y, g, have interpretations. The coordinate 
system thus-described is determined when we fix the direction of 

Z 2 • 

The aerodynamic force perpendicular to the projectile axis is, 

F 

*= < ~ J N + iv J 2 ^) M7— (v JJ-J ^) mUr\ 


and in the non-rotating system the lateral force is this quantity 
multiplied by 
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p 

ijvdp 

0 


This equation then becomes 


*^ (1) = (~ J N + ivJ F )m u2 ~k + (vJ XF + lJ s' )mu2 ^ 


(A.31] 


The component of drag in the plane perpendicular to the 

projectile is due to an aerodynamic force acting on the projectile 

( 1 ) 

with the complex number representation, & , together with that 

component of the axial drag which is in the ' Plane. This 

2 772 

component is kJ n U — . The aerodynamic force perpendicular to the 

u d 

projectile has the complex number Representation described 
by: [ ( J D ~ + ivJ F )k + (vJxF + 1 JgH) 2 • 'll 16 rea l axis is 

perpendicular to the trajectory pointing downward and the imaginary 
axis is horizontal pointing to the left. The quantities k and (i 
are the cross velocities divided by the axial velocity U and the 
cross angular velocity multiplied by d/U. 

The final form of the equation is then obtained. For the sake 
of abbreviation, the term p 

fvdp = <P • Making the same transforma¬ 
tions for £;, f), T : 0 

1=1'^. = Uc7e" i<p ) 'P = [k 'U+\ (-J n -gdU~ 2 sinQ) u-iyku] e~ 1( V^ 


= [k '-ivk-k (J n +gdU sin0) ] ——e 


cr -i<p 


(A.32] 


Similarly, 


f| = [|j/-iYH-|i( J D +gdU~ 2 sinQ) ]—-e 2<p 

a 


(A.33) 


and gT = (y'-ivy)-^e 2< ^. 


(A.34) 


These are now inserted into equations A.23 leading to the following 
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equations: 


m[ \ "-A ( J D +gdU 2 sin0) -i^] = -mgT) -^-e^ < * > =C 1 A — +C 2 |i— -iffy— 5 - 


U U '2 


2 d d 

B[\i'-\i(J D +gdU sin0) ] -iAv^=c 3 A —+C 4 ix— 


(A.35) 


Y'=ign, g ± '= (yil-yn) — 


■B 2 

Substituting in the constants C 4 . . . C 4 , G lt -- =k (the transverse 

radius of gyration in calibers) gives: m< ^ 

A = { JJ lv JgdU sin0) A + g+i) \x~y du 


_o —2 —2 —2 V 

\i'= {-vJ T ~iJ M ) \k +(J D ~k J H +ivk J XT +gdU sin0+iA—)|i 


Y'=isrH/ g 1 '=(YJi-YJA) 


(A.36; 


Finally, the solution of the equation of motion will be 
completed by using some approximations which will simplify 
equations A.36: 

1) Each J-term contains a density factor that is very small in 

-4 

magnitude, on the order of 10 

2) A/B will be approximately 1/10. 

ad -3 

3) As with the density factor, -^-<10 

4) Derivatives of the J-terms ar^ negligible. 

The equations A.36 will be abbreviated as follows: 

A' = a 1 A+a 2 H + jb 

I x' = a 3 A+a 4 jA (A.37) 

with the constants equal to: 

a l =J D~ J N +1VJ E 
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a 2 ~ VJ XF +1 

( — V J ■fc 

a 4~ (J D ~k~ 2 J H +gdU 2 sin0+i^v] 
i>= (gXsin0-y) du 2 . 


Next, we will eliminate H from equations A.37 by substituting 
the first into the second equation: 


A."+X'(-a 1 -a 4 - 


l 2 a 2 
— ') +A. (a-]_a 4 -a 2 a 3 +a-]_— 


a i ,+i< ^r +a 4'i 


To deal with the derivative term a 2 ', recall equation A.30. The 
a 2 ' 2 

terms — , a 1 ' are of the order of J , whereas the other terms in 
the coefficients are of the order of J. So, they can be neglected. 
The result is: 


X "+X ' (-a 1 -a 4 ) +X {a 1 a A -a 2 a- 2 ') + b(a 4 --r ') =0 


(A.38) 


The solution of the homogeneous part of A.38 is done with a 
change of variable to eliminate the first derivative term (the 
particular solution will be addressed later): 


\f (a 1 +a 4 )dp 


X - ge 


(A.39) 


From this equation, the two derivatives X '= [q'+-^ (a 4 +a 4 ) g] e a and 
A. "[g"+(a x +a 4 ) g"+^-(a-L+a,^) 2 g+^-(a 1 +a 4 ) 'g] e a are found. 
Substituting these into A.38 gives: 


q"-z 2 q = 0 


(A.40: 


where r 

= — [ (a 1 +a 4 ) 2 -4 (a 1 a 4 -a 2 a 3 ) ~2 (a-j_ + a 4 ) ] =— [ (a]_~a 4 ) +4a 2 a 3 ~2 (a-^+a 4 ) ] 


72 



To find the approximate solution of equation A.40, let 

Zs(lng) '=—g'; Z' = -^-q' 2 + — '' = -Z 2 +i 2 . Then, Z'+Z 2 -r 2 =0. If r 
Of q o 

were constant then tne equation would have a solution Z=±r. r is 

not really constant because of the v in the a 1 ...a 4 terms, but r 

is nearly constant which allows a trial solution of Z=r+e,e<r. It 

2 X ' 2 

follows e "+2re+e +x'=Q. Let £=-—— then neglecting s', e 


quantities satisfies the equation approximately. The error is of 
2 2 

the order {s'+s )/r . The solution is then: 


Z~r~- 


and hence, 




Using the other sign in Z=±r, 


/t-s? 


(A.41! 


(A.42) 


g = e 


Solving for A, we also expect two solutions 


j[a 1 +a 4 -— )/2±r]dp 


(A.43) 


A = e 


Now, to express the coefficients in A.43 in terms of aerodynamic 

coefficients, the coefficients ^2. +a 4 ^ rom A.37 can be summed and 
. 2 

also substituted into the value of 4r found in A.40. The result 
is aided by using A.30 and yields: 


4 X - A -+4k J M + [J N Jjj 

b a 


-k 2 Jrr-{2Jr r ~J 


. . md‘ 


2 iAV 


(A.44) 


Several specialty terms are required to complete the 
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computation of 


The stability factor, 


A 2 v 2 


4BjC Jr 


0 N 1_ 5 


as well as 


(A.45] 


„ ,222 

R(4r 2 )= A . 

S 2 
r ' 

Neglecting the imaginary part, if — is computed precisely, then 

the denominator multiplied by its conjugate complex with a 

magnitude comparison leads to the same result as if neglecting the 

o 

imaginary part of 4r (except for r«l which is a special case not 
discussed). 


_'= (4r2> (vV) '-.n+n ,j d -j a b*1+ 3 d u- 2 sin0*-h. 

r 2 (4r 2 ) 2 v 2 o 2 v ° ° A A 


(A.46; 


The yawing motion is then a combination of the two solutions 


o ^[J B -J N -*'V J A £ r +iv 4 )±( -p v2t4k ' 2 V( J »- J D-*'V( 2J r- J A>-^) i T 5) 2 1 d P (A. 47) 


2 e “ 


An approximation for the solution of a stable projectile is 

'2 -j^ 

such that the ratio of the imaginary and real parts of i—<— for 

£ U 

stable projectiles. Using the binomial theorem to express the 
square root, 


‘=a+b - 2r=i — va+ [J N ~J D ~k z J H ~ (2J T ~J A )—— ]/a. 

ID /i. 


(A.48: 


In that case, near stability the yaw motion is given by a 
combination of the two solutions: 
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o •§/ [ (J D - J »-*' 2j H +J i. i!! i- )±<J »- J D - ic ' 2 - (2J r- J >i> >/° + nr <1±0) 1 dp 



(a.49; 


This was the form that was used at Aberdeen Proving Grounds for 
reduction of data. An alternate form for a stable projectile can 
be written as follows: 


/ , . 2. / ^ 2 / 222 i2 

r (4r ) A vv a v a . , md , 2 

— =-2~“--T- ( "-5 — ) (J D~ J A—r^ a ' 

1 2(4r ) B Z B Z A 


(A.50: 


and then using equations A.43, A.44, and A.48 leads to: (A.51) 

■1 j [2 Jjj-Jjj-Jc -2 J j/ +gdcr" 2 sin0- (J D -md 2 ~ +gdC7" 2 9in0) /o 2 ± J T -J A ) ] /a* [ivA (l±o) B] ] dp 


The criteria for stability can now be adequately discussed. A 
projectile is deemed to be stable if small disturbances have no 
permanent effect, that is, if the yawing motion does not, in the 
limit, depend on the initial yaw and the initial yawing motion. 
The yawing motion of a projectile is given as a linear combination 
of the two solutions of A.43 plus a particular solution, Xp so that 
X “1,^ X 2 . ^ X, 2 X ^X p . 

The projectile is stable if X 1 , A 2 decrease in magnitude so 
that in the limit only the particular solution Xp remains (it 
becomes independent of the other terms). The particular solution 
need not be-small though, for a large spin, the rate of precession 
is small, and the direction of its axis will change slowly. Since 
the trajectory curves downward after the apex of flight, the yaw 
will become larger. The question then becomes, "what are the 
necessary and sufficient conditions that the solutions given by 
A.47 approach zero as the arc length, p, increases?" One can see 
that the yaw will decrease if and only if the real part of the 
first bracket under the integral is negative and in absolute value 
greater than the real part of the second. 
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Let 


a-( J D +J N +k j h ' 

A 2 V 2 -2 
b=~— ^ +4ic Z J M , 


(A. 52 : 


A ”2 /nd 

c=2 — v [J N -J D -k J H - (2J t -J a ) —— 1 SO that 


q 2 J 
0 


j" [ (-a+iv^) ± (jb+ic) 2 ] dp 


The condition of stability is: 


a> IMy/b+icI . 


2 2 

Reformulating the stability criterion b+ic=\Jb^+c (cos<J)+isin(j)) 


where cosd)= — ^ and sin(t>= ——— 


From De Moivre's theorem: 


^ 2 + c 2 


b 2 + c 2 


j+lcl =4^b 2 +c 2 -|( 1 + -=-) =* (b+\/^ 2+ c 2 ) ■ 

\ 2 N 2 


(A.53: 


So that the stability criteria become: 


a>* (b+\/b 2 +c 2 ) 


where 


-a=R(a 1 +a 4 ) 


b=R4r' 


c=Jm 4r . 


Additional steps are taken in reformulating the stability 
criteria in terms of aerodynamic coefficients: 
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and finally. 


> -| {b+sjb 2 +c 2 ) , 

a>0 

2 a 2 > b+/b 2 +c 2 , 

a> 0 

2a 2 -b > \lb 2 +c 2 , 

a> 0 

{2a 2 -b) 2 > b 2 + c 2 

, a 

4 2.2 

4a -4a b > c , 

a>0 

2 2 4 

-4a b > c -4a , 

a>0 

2 c 2 

~b > -a 2 + , 

4 a 

a>0. 


(A.54) 


2 2 c 2 2 

Noting that a is of the order of J and -a is of the order 

4 

,2 _ t 2 2 ,.2 

A 2 u V A 2 

-v - = - > J, Jb ~ —-v > J, therefore, the a's can be 

B 2 J 2 100 B 2 


neglected. Let 


2 2 
A 2 -2 A -2 mr? A 

b = v +4Jc J M , and c = 2 — v [ J N -J D ~k J H ~ (2 J T ~J A ) —] =2 — v-f. 


B 


A 


B 


such that now we can rewrite inequality A.54 as 
2 _ _ (f24v ) 2 _ 4 B 2 k~ 2 J„ j=2 


A "..2 -2 T . < f2 | V)2 .2 . , f . 

2 V AT 4 a _, 2..2 2 (A.55) 


B 


aAv 


1 4B 2 k~ 2 J M 

Recalling the stability factor, S, — =- allows the 

■* ^ .. ^ 


following steps: 


A v 


/ -r -r , -2 _ T md , _ _ 
a - ( J D +J N +k J H J A~~ S 1 


-2 


md 


a + f—{ J^+J^j+k J u J a—-—+ (Jj\j J p, k J (2Jrp J %) 


md‘ 


H A 


N U D 


/ jj \^U rj-t U 
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a-f-(~J d +J N +k Jh~ J A~~X~ ( J N J D k J H ( 2J T J A> —~ s 2 

These terms S-j_, S 2 , S 3 are very important to the rest of the 
discussion of stability. The final criteria for stability is: 


1 < 

s si 


S 1 > 0 


(A.561 


. -2 _ md‘ 

where S 1 = J N +k J H -J D ~Ja~~X 


5 2 = 2J N 2J d ~2J t , 

5 3 = 2k~ 2 J H +{2J T ~2J A )^~ 


(A.57) 


Using the relationships between the S terms and f, the 
stability criteria can be rewritten. 


S 2 S 3 (S-L + f) (S x -f) s 2 _ f 2 


s 1 >0 


(A.58) 


-|/[ (-S 1 +i2ic 1 Jjrf3)±J5+Ic]dp 


Then, A«e 


-2 , A .2 ,, -2 


where b=4ic J M -——v =4k J M (1- 


a 2 v 2 


, M —v — _ 2 


) = 4k 


4B k 


C=2^Vf = 2^4k 2 J m JS-f 


and because S>0, it follows that S>1 from A. 58. This is the 

"classical stability condition." This is not sufficient though. 

The general motion for .r , 

-i I v dp 

l = Ule J , 


with the two solutions for A represents the superposition of the 
two cyclic motions, each decreasing in amplitude. This is the 
epicyclic motion caused by the coning and spinning simultaneously. 
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For S values around 1.5, the epicyclic yawing motion will have many 
maxima. The successive values of the yaw at maxima may increase 
without bound. 

If the classical condition for stability is not met (S<1) , 
then the final step of the inequality A.58 can not be satisfied. 
The yaw will grow steadily and the motion is similar to the falling 
motion of a top as the spin becomes too small. 

It would be difficult to meet the condition S^CO since the 
terms J N , 10 J D . If J H < 0, then this condition could be met, 
but this is unlikely to happen since it would require a shift of 
the center of mass, which for our purposes is not practical. 

Finally, the stability found in the product of the S 2 , S 3 
terms may be sufficient to cause instability. Recall that: 


1 F F °2 f 

— < l-— for S-, >0, —- could become > l. —— =1 + -—, 
S 2 1 s 2 Sl Sl 



This would happen if either of the S 2 , S 3 terms become negative. 
But, because S 2 + S2 =: 2S 1 , S 2 +S 3 > 0 , always. Depending on the sign of 
the Magnus moment coefficient, Jj ,, a change in magnitude could 
cause either S term to become negative and result in an unstable 
projectile. 
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APPENDIX B 


The laser melting of a projectile surface is the alternative 
to multiple pulse dimples. To heat the surface to a depth of 
greater than 3 mm, the surface temperature profile will appear as 
in Figure 29 below. __ _ _ 



The melting temperature at a depth for steel is related to the 
surface temperature by: 


Tw = —- = 1450 C, T„ = 3900 C. 


(B.l) 


The heat diffusion thickness is: 


e = s.lKtT), where < steel = 0.05——, t D 


(Vi 

2 K 


= 0.45 sec, 







The required heat is best deposited in a time, tf < t D . The 
average temperature to which the material must be raised is 
estimated to be: 


- '^M +r ^s 

T = = 2700 C. 


(B. 2) 


This gives a required fluence of: 

F 0 t f = pe[L M + (T-T q ) C s = 4700 J/Cir> 2 . ( B .3) 

For continuous laser illumination with a sufficiently large 
burn spot (radius 30 cm) on a projectile rotating with a period of 
4 milliseconds, an equal distribution of the fluence is achieved. 

The absorptance of steel is approximately 30% (Schriempf, 
1974) . This means that for a laser intensity, I L , the fluence 
becomes: 

F o t f =tf-0.3-I L = 4.7KJ/cm 2 . (B.4) 


This gives a required laser energy of 


nr‘ 


„ 4700 _ 2 

= 2 x -.3 J/cm , 


(B. 5) 


where the factor 2 takes into account the energy required on both 
sides of the projectile and the factor 0.3 accounts for the 
reflectance. The result of solving equation B.5 for the energy is 
nearly 90 mega joules required. For a tf < 0.45 sec, this 
corresponds to a laser power of 200 megawatts. 

If the burn spot is reduced to an 8 centimeter radius and it 
is assumed the laser energy is moved over the projectile body, then 
substituting values into equation B.5 results in an energy 
requirement of only 6 megajoules and a power of 14 megawatts. 
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